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Abstract 

A method for designing an automatic flight controller for 
short and vertical take off aircraft is presently being developed 
at NASA Ames Research Center. This technique involves transforma- 
tions of nonlinear systems to controllable linear systems and takes 
into account the nonlinearities of the aircraft. In general, the 
transformations cannot always be given in closed form. Using 
partial differential equations, an approximate linear system 
c<illed the modified tangent model, was recently introduced. A 
linear transformation of this tangent model to Brunovsky canonical 
form can be constructed, and from this the linear part (about 
a state space point x^) of an exact transformation for the nonlinear 
system can be found. Here we show that a canonical expansion in 
Lie brackets about the point x^ yields the same modified tangent 


model . 
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I . Introduction 

Suppose we have a nonlinear plant for which we are to design 
a control scheme. For example, George Meyer [1] , [2] , [3] , [4] , [5] , 

[6] at NASA Ames Research Center is presently developing an auto- 
matic flight controller for the UH-lH helicopter which takes into 
account the nonlinearities involved. His design technique depends 
on a theory giving necessary and sufficient conditions for a non- 
linear system to be transformed to a controilable linear system 

[7] , [8] , [9] , [10] , [11] . In other words state and control coordinate 
changes can be implemented to simplify the problem. Thus the 
method is to move 
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from the nonlinear model of the plant to the linear system in 
Brunovsky tl2) canonical form. To have the helicopter fly a 
prescribed trajectory, no gain scheduling is needed in Meyer's 
technique, because through on-line computations of the transformation 
and its inverse, we always see the same trivial linear system. 

Meyer's nonlinear system is in block triangular form, and 
it is not difficult to find a transformation and its inverse. 

However, the transformation theory in [8] applies to systems which 
are much more general than block triangular. Formally, the 
desired transformations can be constructed by considering a system 
of partial differential equations, which can be reduced to ordinary 
differential equations. It is not always possible to solve such 
equations in closed form, but cases where this can be accomplished 
are presented in the Ph.D. thesis of H. Ford 1131. Numerical 
techniques for constructing approximate transformations in certain 
situations (e.g. under the conditions due to Brockett [14]) are 
introduced in [131. 

It is appropriate to develop a method to build an approximate 
transformation in all cases. In [15] wo considered this problem 
in view of the partial differential equations from [8]. We found 
a related set of partial differential equations, the solution of 
which yielded a linear approximation to an actual transformation. 

In fact, given a point Xq in state space, we are able to construct 
a mapping which is the linear part of an actua’ transformation about Xq , 
without knowing the transformation itself. In this process we 
introduce an approximating linear system called the modified tangent 
model. If one is working around an equilibruim point of the drift 
term in our nonlinear system, this modified tangent model is the 
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linear system we obtain by truncating the Taylor series of our 
nonlinear system about the point. However, if we are away from such 
an equilibruim point, the modified tangent model can provide a 
different system than the usual linearisation (an example is 
presented in [15]). If one is interested in tracking a certain 
trajectory, then modified tangent models are constructed at various 
points along the trajectory. 

The purpose of this paper is to show how one obtains the 
modified tangent model by an expansion which parallels, but is 
unequal to in most cases, the Taylor approach. This technique 
is included in an article [16] by the authors giving canonical 
forms and canonical expansions for nonlinear systems. 

We show that the partial differential equations and canonical 
expansion yield the same approximate controllable linear system. 
Moreover, the form of the approximate linear system and knowledge 
about nonlinear control theory convince us that it is the appropriate 
model for designing a controller for many nonlinear plants. 

II. Expansions and the Tangent Model 

In order that a nonlinear system be transformable to a linear 
system as in [8] , it must be of the form 

m 

(1) x(t.) = f(x(t)) + u i (t)g i (x(t) ) , 

where f , g g are C vector fields on IR n , and g.,g_, . . . ,g m 
are linearly independent (this is assumed for convenience) . Now 
the results from [8] are local in nature (global theorems are 
presented in [9]), but for the sake of simplicity we assume that 
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the system is transformable to a controllal$¥e linear system on 
all of R n . This linear system is in Brunovsky form, 
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(2) £(t) * A Q y + B q v. 

and with Kronecker indices k. > k_ > . . 

l — i — — m 

For vector fields f and g, (f,g] denotes the well known 
Lie bracket, and 


(ad°f,g) = g 
(ad 1 f,g) = [f,g] 
(ad 2 f,g) = [f, [f ,g] ] 


(ad k f,g) = [f, (ad k_1 f,g) ] . 

We define 

C = { , [ f , g A 1 » ... , (ad' 1 f ,gj^) ,g 2 , (f ,g 2 ] , ... , (ad 2 f,g 2 ) 

... ,g m ,[f.9 m j, ... , (ad K m” 1 f,g in ) } 

K -2 K -2 

Cj ={g 1 ,[f,g 1 l, ... , (ad j f ,g x ) ,g 2 , [f ,g 2 ) , ... ,(ad'j f ,g 2 ) , 

... ,9^ » »9 In 1/ ... > (ad j f , 9 ^ )}for j = l,2,...,m 
In [8] it is shown that system (1) is transformable to system (2) 
if and only if (with possible reordering of gj.,g 2 , ... ,g m ) 

1) the n vector fields inC are 1 inear Jy independent, 

2) the sets Cj are involutive for j=l , 2 , . . . ,m, and 

3) the span of equals the span of C. (1 C for j=l,2,...,m 


We assume that our system (1) satisfies these three condi- 



tions. Using the partial differential equations approach from 
[8] we introduced the modified tangent model about a point x^ in 
x-space 

(3) x(t) = f(Xg) “ AXq + Ax + Bu. 

It is shown in [15] how to construct an approximate transformation 
using this model. Here A is an nxn matrix and B = (b^,b 2 # ••• ,b m ) 
is an m tuple of n vectors that satisfy the equations (take + 
for k even and - for k odd) 

A 1 ^ = + (ad k f (x Q ) , k = 0,1, — ,k ^ 

(4) A k b 2 f + (ad k f,g 2 ) (x Q ) , k = 0,1, ... ,< 2 

A k j = + (ad k f ,g ro ) (x Q ) , k = 0,1, ... , 1 ^. 
m 

Equations (4) are nonlinear, but there is a simple method for 
computing A and B. Let D be the set of vector fields 
{ (ad K ^-f ,g 1 ) (x Q ) ,ad Kl “ 1 f ,g ][ ) (x Q ) , ... , (ad <2 £ ,q^) (x Q ) , 

(ad^ 2 f ,g 2 ) ( x Q ) , (ad K2 “ 1 f ,g 1 ) (x Q ) , (ad K2 ~ l f ,g 2 ) (Xq) , ... , (ad K3 f ,g 3 ) (x Q ) , 
... ,g (x Q ), g 2 (x Q ), ... ,g m (x 0 )}. Before forming this set 
checks such as = < 2 or > k^, etc, should be made and no 
duplications should be included. 

Now we introduce an interesting (n+m)x(n+m) matrix E. Let 

Yl 1 

the first column be (ad i f,g^) (x^) followed by m zeroes, the second 

column will be the second element of D followed by m zeroes, ... , 

the n— column be the n— element of D followed by m zeros, the 
til 

(n+1)-— column be g^(Xg) and m zeros, ... , the last column be 

g (x rt ) and m zeros. 

^m 0 

Ignoring the last m components, the first column of E is + A^b^, 
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the second + A K 1 ^b , , ... , the n— + Ab # ... . and the last b . 

— — m m 

is shown in [ 13 ] that there is an orthogonal coordinate change on R* 
so that all entries above the first m superdiagonals are zero 
and the elements below the n— row remain unchanged. For our 
to ke purpose we can assume that E is initally in this "generalized lower 
St»i4 ke** Hessenberg" form, because knowing A and b in these coordinate, 

we can return to the original A and b through an orthogonal change 
of coordinates. 

Hence we have by (4) 



where * indicates a possible nonzero entry (recall that 
9 m' g m-l' '9 1 were assumed to be linearly independent) and 

the first * in b^ is in the (n-m + 1)— row. 


We examine from (4) 



k 
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with the first * in the column being at the (n-m) — level. This is 


a a. _ 

11 12 

a 21 a 22 




and we easily compute a i n * a 2n * ••• • a nn * Similarly, 

^m-l “ ~ tf ^m-l 1 (x 0 ) yields a l (n-1) ,a 2 (n-1) ' *** ,a n(n-l) * 

Abj. = -|f.g 1 l(x 0 ) yxeids a 1(n _ m+J) f a 2 (n_ m+ i y ••• ' a n ( n -m+l) * 

Next we consider 


» * * # 


A 2 b * (ad 2 f,q (x n ) 
iu m u 

if the vector field on the right is in the set D. Writing this 
as 


a\ = A(M, ra > - («d 2 f„ 1ra )(x 0 ), 

*> 

and knowing Ab, and (ad‘f,q ) (x rt ) , we can compute a,, . ^a*, „v , 

ro mu i(n~nw 

... » a n ( n . m j • Continuing in this way we can solve for every entry 
in A, and the method of solution is readily implemented on a 
computer (or by hand) . We remark that the above equations can be 
solved because of our assumption that, the set of vector fields in 
the set C are linearly independent. □ 

Now we show the canonical expansions like those in (161 give 
us the modified tangent model. First we rewrite the set C so that 
the vector fields appear in a different order. 
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C * { (ad K l 1 f,g 1 ) , (ad Kl ” 2 f ( g x ) , ... , (ad* 2 " 1 ! ,g x ) , (ad K 2” 1 f ,g 2 ) 

(ad* 2 ~ 2 f #gi ) , (ad* 2 “ 2 f ,g 2 ) , . . . , (ad* 3 " 1 ! ,g x ) , (ad K ‘ 3 “ 1 f ,g 2 > , (ad* 3 ’ 1 ! , g 3 ) 


(ad ^ f#g^)#(ad 3 f,g 2 ) f (ad K: 3 f,g 3 ), ... ,g 1 # g 2 , ... # g ro >. 

Note that Kj > or K ^ E etc » should be checked before this 

set is formed. For example if * k- 2 , the first element is 

(ad k l 1 f,g^) , the second (ad K 2 ’*f ,g 2 ) , etc. 

We introduce new independent variables s.,s,, f ... ,s such 

12 n 

that Sq * ^ 3 l 0 ' s 20 ' ... »s n Q> is the same point as 

X 0 = ^ X 10' X 20' *•* ' x n0^ ls our state space. The parameter is 
along the integral curve 3 of the first vector field (ad* 1 - 1 f , 9 ^) 

in C , Sj is along the integral curves of the second vector field 

in C, s 3 is along the curves of the third, ... , and s n is along the 

integral curves of g m - Notice that our system (1) is linear 

in the controls u 3 ,u 2 » ... ,u m and these control variables are 

treated at the same level at the s.,s_, ... ,s . That is. in our 

1 * n 

linearization we do not want terms like u^s^,u^s 2 , etc, because 
these are not considered to be linear. In the following process 
all terms of degree greater than one in u.,u_, ... ,u ,s.,s.,, ... ,s 
are included in the notation + ... . In taking infinite expansions 
f and g are required to be real analytic, but we are only interested 
in finite truncations of such expansions. 

First we rewrite our nonlinear system (1) 


(51 x - * s 2 ' ... , s^) + 92 ^ s i ,S 2' • • • » s^) u 2 + ... + 

^ s i ' s ? ' *•' ' u m + f(s.,s_, ... ,s ) . 


Now we expand in the Sj variable along the integral curves of 
1-1 


(ad 


f,d,> 
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* " ^1 ^*10 * a 2 ' * * * ' 8 n^ U 1 + g 2 ( *10 #S 2 ' * * * ,8 n^ u 2 + * * ’ + 
g m* s l0' s 20' *” ,8 n^ u m + f * s i0' 8 2' *** ,8 n^ 4 
[f# (ad f #g^) ] ^ s io» s 2‘ ••• * ^ 8 l” 8 10^ ^ *•• 

Expansion along the integral curves of the second vector field in 
C, which we assume is (ad 1 “ f,gj) , gives 

x - t, j^ s ^o ,s 20' *" ' 8 n^ U l + g 2^ s 10' s 20' **' ' s n^ u 2 + • • • + 

g m (s 10' s 20 s n )u m + f(s 10' s 20 S n ) + 

n 2 

( f t ( ad f 1 ^ ^10 1 ^20 * ••• i s^) (82"®2 q) ^ 

(f , (ad K 1-1 f .g^ ] ( S 10 ' S 20' * * * ' S n^ ^ S l -S 10^ + 

Continuinq in this way, we arrive at our last step which is an 
expansion in the s n variable that provides ( with s Q *= 

* s 10' s 20' * * * ,8 n0 ) ) 

X » g 1 (s 0 )u 1 + g 2 (s 0 )u 2 + ... + g m (s 0 )u m + f(s 0 ) + 

(6) (ad Kj -f ,gj ) (s Q ) (s^-Sj^q) + (ad*' 1 ' 1 f , g^) ( s Q ) ( s 2 -s 2Q ) + ... + 

lf ' g m Us 0 )( V S n0 ) + ••• 

Since Sq corresponds to Xq, the important Lie brackets are 
g l (x o) ' g 2 (x 0 J ' *•' ' g m ^ x 0 ^ ^ f ' g 1 1 ( x o ' ' 1 f * 9 2 ) ( Xq ) » ••• » (ad K (x Q ) . 

Another way to find this set is to take the elements of C evaluated 
at Xq plus (ad m f ,g m ) (x 0 ) , (ad^^^f ,g m _ 1 ) (x Q ) , . . . , (ad K ,9^) (x 0 ) . 
Thus if we wish to find a linear system 

( 7 ) x * f(Xjj) - AXq + Ax +Bu 

that emphasizes the linear part of the system in (6) we would need 
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to solve 

A c + ( sd f»9^) (^g) » ^ B 0,1, • • • ^ ^ 
“ + (ad k f ,g 2 ) (x Q ) , k = 0,1, ... ,< 2 


A k b m = + (ad k f ,g m ) (x Q ) , k = 0,1, ... ,< m ) 
where B = (b^.b^ ... ,b m ) and + is for k even and - is for k odd. 

Thus (7) is exactly the modified tangent model that we defined 
in terms of the partial differential equations. 

III. Conclusion 

Given a nonlinear system (1) we can use expansions of the 
system about a point Xq in terms of variables associated with the 
vector fields in the set C to produce an approximate linear system. 
This resulting linear system is the modified tangent model that 
is found by considering the partial differential equations that are 
solved in constructing an exact transformation of the nonlinear 
system to a controllable linear system. 
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